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Abstract 

One of the useful and practical methods for solving quantum-mechanical many-body systems is 
to recast the full problem into a form of the effective interaction acting within a model space of 
2 ■ tractable size. Many of the effective-interaction theories in nuclear physics have been formulated 
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by use of the so called Q box introduced by Kuo et.al. It has been one of the central problems how 



J> ■ to calculate the Q box accurately and efficiently. We first show that, introducing new basis states, 
00" 

CN ■ the Hamiltonian is transformed to a block-tridiagonal form in terms of submatrices with small 

o. 

. dimension. With this transformed Hamiltonian, we next prove that the Q box can be expressed 



in two ways: One is a form of continued fraction and the other is a simple series expansion up 
to second order with respect to renormalized vertices and propagators. This procedure ensures to 
derive an exact Q box, if the calculation converges as the dimension of the Hilbert space tends to 



H . infinity. The Q box given in this study corresponds to a non-perturbative solution for the energy- 

dependent effective interaction which is often referred to as the Bloch-Horowitz or the Feshbach 
form. By applying the Z-box approach based on the Q box proposed previously, we introduce a 
graphical method for solving the eigenvalue problem of the Hamiltonian. The present approach 
has a possibility of resolving many of the difficulties encountered in the effective-interaction theory. 

PACS numbers: 21.60.De, 24.10. Cn, 02.30.Vv, 02.60. Cb 
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I. INTRODUCTION 



In nuclear many-body physics various methods have been proposed to solve the 
Schroedinger equations for nuclear many-body systems starting with realistic nucleon- 



nucleon interactions. These methods, which are ca 
;he Green's function Monte Carlo (GFMC) method 



3,3, 



led the ab initio calculations, include 
2j, the no-core shell-model (NCSM) 



the effective interaction for hyperspherical harmonics (EIHH) method j^, the cou- 
pled cluster method (COM) 6-8| and the unitary-model-operator approach (UMOA) js-ll]. 
Much effort has been made also to diagonalize a matrix of a many-body shell-model Hamil- 
tonian in a huge dimensional Hilbert space on the basis ofQ, 13], or alternatively to[l^ the 
Lanczos method. 

The effective-interaction theory aims at recasting a full many-body problem into the 
form of a small-dimensional model-space problem with a corresponding effective interaction. 
Many studies on this subject have been made both as regards formal theories and their 
-2l|. This effective-interaction method has been developed to apply to new 



applications 
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22, 



23j and many-boson systems 24 1. 



fields of many-body physics such as quantum dots 

Most of the effective-interaction theories given to date have been formulated in terms 
of the Q box introduced by Kuo and his collaborators 
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. Originally the Q box 

has been defined as the sum of linked and unfolded diagrams [27|]. In the algebraic or 
non-diagrammatical approach the Q box is equivalent to the energy- dependent effective 



interaction given by Bloch and Horowitz 28 



and Feshbach 



30 



29[ which has been studied 



extensively on the Brillouin-Wigner perturbation theory 

It has been established that the effective interaction can be expressed as a series expansion 
in terms of the Q box and its energy derivatives. The series can be summed up by using either 
the Krenciglowa-Kuo (KK) [32] or the Lee-Suzuki(LS) 33l-l35| methods. It has been known 
that, in general, two methods have different convergence properties: Many of the numerical 
calculations have shown that the KK method yields the eigenvalues for the eigenstates which 
have the largest overlaps with the chosen model space. However, it has been pointed out 
that the rigorous convergence condition for the KK method has not yet been clarified js^]. 
On the other hand the LS method reproduces the eigenvalues which lie closest to the chosen 
unperturbed energy. Both of the two approaches reproduce only certain of the eigenvalues 
of the original Hamiltonian. This restriction is not, in general, desirable. 
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Another difficulty encountered in actual calculations is the pole problem. The Q box 
itself has poles at the energies which are the eigenvalues of QHQ, where Q is the projection 
operator onto the complement {Q space) of the model space (P space). The presence of the 

poles causes often instability in numerical calculations. Three of the present authors and 

i 

Fujii [38] have shown that it was indeed possible to resolve these difficulties by introducing 
a new vertex function Z(E), called the Z box. The Z-hox approach based on the Q box 
and its derivatives may have a possibility of resolving many of the difficulties encountered 
in the effective-interaction theory. 

At present the most important remaining task would be to establish a method of how 
to calculate the Q box rigorously and efficiently. The perturbative calculation method for 



13, 



27 



32l |. but a basic problem of the 



the Q box has been established and applied widely 
convergence of its perturbation expansion has not been made clear for general cases. Main 
concern of the present article is to propose a non-perturbative method for calculating the Q 
box. 

The formulation in the present study consists mainly of two parts: The ffist one is to 
transform the Hamiltonian to a block-tridiagonal form, where the dimensions of the block 
submatrices are taken to be equal-to-or-less-than the dimension of the P space. With the 
block-tridiagonalized Hamiltonian, the next step is to derive a set of coupled equations for 
determining the Q box. We show that the coupled equations can be solved by employing 
two different recursion methods: The ffist solution is represented in a form of continued 
fraction, and the second one is expressed as a sum of terms up to second order with respect 
to renormalized vertices and propagators. In both of the methods the calculation of the Q 
box can be carried out without matrix inversion of QHQ which is usually a huge-dimensional 
matrix. All the procedures for obtaining the Q box are reduced to calculations of small- 
dimensional submatrices in the block-tridiagonalized Hamiltonian. 

Regarding the block tridiagonalization of the Hamiltonian, the present approach has a 
common aspect to the so called block Lanczos method based on the theory of the Krylov 
subspaces ISj. For a given model space P and a Hamiltonian H, the subspaces leading 
to a block-tridiagonal form of H are determined uniquely. Therefore, the subspaces given 
in the present study are the same as those of Krylov. However, the choice of basis states 
of each subspace is ambiguous. For determining the basis states we employ a different 
calculation procedure from the usual one in the block Lanczos method. Different basis 
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states are introduced, and we show that they are suitable for the purpose of calculating not 
only the Q box but also the eigenstates of H. 

The construction of the present article is as follows: In Section II some basic elements of 
the effective-interaction theory are reviewed. Section III is devoted to rigorous calculation 
of the Q box: In the first step, the Hamiltonian is transformed to a block-tridiagonal form. 
Next we derive a set of coupled equations for determining the Q box. The equations are 
solved to derive two forms of the solutions for the Q box, where one is represented in a form of 
continued fraction and the other is expressed as a sum of terms with renormalized vertices 
and propagators. In Section IV, a method is given for the problem of how to calculate 
eigenstates of H within the framework of the effective-interaction theory. In Section V a 
short review of the Z-hox theory is given. In Section VI, by applying the Z-hox theory, we 
make a numerical calculation with a model Hamiltonian to assess the present approach. We 
propose a graphical method and show that it works well for finding the eigenvalues of H. 
Summary of the present study and some remarks are given in the last Section. In Appendices 
A and B the derivatives of the Q box are given, which are necessary for calculating the Z 
box. 



II. EFFECTIVE-INTERACTION THEORY BY MEANS OF SIMILARITY 
TRANSFORMATION 

Let us begin with a Hamiltonian H defined in a Hilbert space. We divide the space into 
a model space (P space) and its complementary space {Q space). When all the eigenvalues 
of an operator H^s given in the P space coincide with those of H, we call H^g an effective 
Hamiltonian. In the following, we do not impose any particular conditions on H and states 
belonging to the P space nor assume degeneracy of their unperturbed energies. 

There are various ways of constructing Hcs- We adopt the following standard one. First 
we introduce an operator u that maps states in the P space and those in the Q space to 
each other, with the properties [s^ . 

uj = QujP, (II.l) 
cj" = (n > 2). (H.2) 
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The operator uj defines a similarity transformation of if, 

H = e-'^He''. (II.3) 

This reduces to 

H = {l-uj)H{l + uj) (II.4) 

by virtue of Eq.(^^. 

The condition that PHP be a model-space effective Hamiltonian ifeff is that H is decou- 
pled between the P and Q spaces as 

QHP = 0. (II.5) 

This condition is rewritten as 

QHP + QHQcu - ujPHP - cuPHQu = 0, (II.6) 



with the aid of 
different way 



qs. flll.ip and (III.4|) . This equation for u was first derived by Okubo in a 



. Once a solution u to Eq.f lII.6|) is given, i^efr is written as 

H^s = PHP + PHQuj. (II. 7) 

Dividing PHP into the unperturbed part PHqP and the interaction PVP, we write 

PHP = PHqP + PVP. (II. 8) 

The model-space effective interaction Ves is defined as 

V,s = H,s - PHoP 

= PVP + PHQuj. (II.9) 

From the definition of i^es and V^g we see that a central part of determining them is to find 
a solution for u in Eq.f lII.6|) . 



Since Eq.( III.6P is a nonlinear matrix equation for w, it is difficult to find a general solution. 
The following formal solution, however, has been known and is enough for applications. We 
rewrite Eq. (lII.6p as 

QHP + QHQuj -uH^s = 0, (11.10) 
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using Eq. (lII.7p . Here the eigenvalue equation for Hes is given by 

H,s\<Pk) = Ek\<Pk). (11.11) 



If the operator a; is a solution to Eq. (lll.6p . we can verify that eigenstates {|0fc)} belong to the 



P space and each eigenvalue coincides with one of those of H. The effective Hamiltonian 
ifgff is not Hermitian in general; the eigenstates {\4>k), k = 1,2, ■ ■ ■ ,d} are not orthogonal 
to each other. Then we introduce the adjoint states {{(pk\, k = 1,2, ■ ■ ■ ,d} according to the 
biorthogonality condition 



) = 5kk', (11.12) 



where d is the dimension of the P space. The projection operator onto the P space is written 
as 



k=l 



Then, using Eqs. (llLT0|) - (llLT3|) . lo is given by 



1 

= E 7TT77^QHP\<Pk){<Pk\, (11.14) 



and from Eq.( III.7p Hes becomes 

1 

H,^ = PHP + Y,PHQ^-—^QHP\ct>u){M- (H.IS) 

Here we introduce an operator in the P space called the Q box 

Q{E) = PHP + PHQ^-^j^QHP, (11.16) 

where E is an energy variable. The Q box thus defined i s eq uivalent to the energy- dependent 



28| and / or the Feshbach 29|] forms 



effective Hamiltonian referred to as the Bloch-Horowitz 
In terms of Q{E), Hgs is expressed as 

d 

H,^ = Y,Q{Ek)\<Pk){4>kV (11.17) 

k=l 

from which the following self- consistent equation can be derived 

Q{Ek)\<Pk)=Ek\<Pk)- (11.18) 
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The i^eff in Eq. (III.17|) is just a formal solution in the sense that unknown E^, \(j)k), and 
appear in the right-hand side, but the following method of solving is available: In order 
that the solutions to Eq. (III. IIP coincide with those given by Eq. (111.171) , they selfconsistently 
satisfy the iterative equation 

where 4"^'^ and 101"+'^) are the {n + l)-th order eigenvalue and eigenstate of the Q box, 
respectively, given by the n-th order eigenvalue Ej^\ There have been a lot of studies about 



the convergence of this iterative method 



32h37|| . But the condition of convergence is rather 



complicated and it has been known that only some specific solutions are obtained. 

In addition, Q{E) has poles at energies {eq}, where Eq is one of the eigenvalues of QHQ, 

QHQ\q)=eq\q). (11.20) 

These singularities of the Q box lead to some difficulties in numerical calculations 
These arguments suggest that some further improvements are desired for the Q-box method 
although it has been applied widely to practical problems. 



III. CALCULATION OF THE Q BOX BY MEANS OF RECURRENCE RELA- 
TIONS 

Most of the effective-interaction theories formulated so far are based on the Q box. The 
Q box has been calculated on the perturbative expansion methods, but their convergence 
properties and accuracies have not been well understood yet. This is because, as a matter 
of fact, it is impossible to solve the eigenvalue problem of QHQ nor to calculate the inverse 
of {E — QHQ) when the dimension of the Q space is huge. The accuracy of the Q box 
determines that of H^^ and Veg, because errors that arise in the calculations of operators 
and/or matrices in the P space with small dimension are considered to be negligible. 

In the following subsections we describe a method of how to calculate accurately and effi- 
ciently the Q box. We first transform if to a block-tridiagonal form. With this transformed 
Hamiltonian we derive a set of coupled equations for determining the operator u. We shall 
show that these coupled equations can be solved in two ways by introducing two types of 
recurrence relations. The properties of two solutions for the Q box are discussed. 
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A. Block tridiagonalization of Hamiltonian 



We transform the Hamiltonian H into a tractable form by changing basis vectors. First 
we introduce 



The Yp is an operator in the P space, which is Hermitian and positive semi-definite, that 
is, y^}^ > in the eigenvalue equation 



Suppose that di eigenvalues are nonzero among {y^ }. In terms of the eigenvectors {\pk), k = 
1, 2, ■ ■ ■ ,di} with nonzero eigenvalues, we define normalized vectors {lo'^^'')} in the Q space 
as 



They are orthogonal to each other and span the (ii-dimensional subspace Qi in the Q space. 
Then the projection operator onto the Qi space becomes 

Qi = El^^)(?i'Y (ni.4) 

k=l 

The complement of the Qi space in the Q space is defined as 



Yp = PHQ ■ QHP. 



(III.l) 



yp\pk) = Vk \pk)- 



(III.2) 




(III.3) 



Q,=Q-Qi. 



(III.5) 



Equation (1III.3I) indicates that 




(III.6) 



k=l 



then we have 



QHP 



QiHP 



(III.7) 



which leads to 



Q^HP = 0. 



(III.8) 
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a= Q-Q^ 



FIG. 1: The image H{P) by the mapping H. 



Thus the image H{P) by the mapping H is given as a sum of the P and Qi spaces as 
depicted in FiglH 

Next, a similar manipulation with replacing P and Q with Qi and Q^, respectively, leads 
to another orthogonal system. We introduce 

Yq,=Q,HQ,-Q,HQu (III.9) 

and write its eigenvalue equation as 

)=2/f )• (III.IO) 

The eigenvectors {l^'^^^)} belong to the Qi space and accordingly are given as linear com- 
binations of ilq^k'')^^ ~ l'^,--- ,di} in Eq. flIII.3p . Suppose also that d2 eigenvalues are 
nonzero among {y^ }. New orthogonal bases 

> = ^^OiffQilft"'), k = 1,2,- ■ -,112 (m.ii) 

are derived. The (i2-dimensional subspace Q2 is defined by them and the projection operator 
onto the Q2 space is expressed as 

d2 



Q2 = $^|gf)(gf|. (III.12) 

k=l 
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The projection operator Q2 has the properties 



Q2HP = 0, (III.13) 
QiZ/Qi = Q2HQ^, (III.14) 

g2^gi = o, (III. 15) 

where Q21 the complementary space to Qi + Q2 in the Q space, is written as 

Q2^Q-Qi- Q2. (III.16) 
Repeating these manipulations leads to the following: Decompose the Q space as 

Q = Qi + Q2 + --- + Qn + --- ■ (111.17) 

Basis vectors of a subspace Qm, namely, {\qk"^),k = 1,2,- ■■ ,djn}, define the projection 
operator 

dm 

Qm-Y.\^tWr\ (ni.i8) 

k=l 

The basis vectors are given as follows: Introduce Yq^_^ as 

Yq^_, - Q^_,HQ^_, ■ Q^_,HQ^^, (III.19) 

with 

Qm-l = <5 - (<5l + <52 + • • • + Qm-l). (111.20) 

Its eigenvalue equation is 

YQm-M't-'^)-y^r^\ci't-'^)- (ni.2i) 

In general new orthogonal bases 

\<lt^) = -^Qm-iHQm-Aq't-') (III.22) 

are derived from the eigenvectors {l?^*"'""^^)} with nonzero eigenvalues {y^^}- They span 
the subspace Qm- When all the eigenvalues {y^^^} are zero, the procedure ends because the 
eigenstates of H reside in the subspace P + Qi + Q2-\ h Qm-i- Here we note that we are 
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not interested in any eigenstates that are decoupled from the states in the P space. With 
the projection operators Qm and Qm-i we obtain, from Eq. pil.22|) . an expression written 
as 



(III.23) 



We conclude from the above discussion that 

PHQm = QmHP = Q (m>2), (111.24) 

QmHQm+k=Qm+kHQm = Q {k > 2) (III.25) 

hold for the subspaces {P, Qi, Q2, ■ ■ ■ , Qm, ■ ■ ■ }• This means that the given Hamiltonian H 
is transformed to a block-tridiagonal matrix 

f PHP PHQi ■ ■ ■ \ 

QiHP QiHQi Q1HQ2 

Q2HQi Q2HQ2 Q2HQS ■■■ 

Q,HQ2 Q^HQ, ■ ■ ■ 

'■■ '■■ ■■ ■■ J 

where each block matrix is at most d-dimensional. Thus the image H{Qm) by the mapping 

if is a sum of adjacent subspaces Qm-i, Qmi and Qm+i as depicted in Figj2j 



H 



\ 



(III.26) 




















p 








0^ 




J 



FIG. 2: The image H{Qm) by the mapping H for m > 2. 

From FigsiH and [2] it is easy to see that the image of the mapping H of the P space 
becomes 



H{P) = P + Qi. 



(IIL27) 
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The image of the successive mapping is given by 



H\P)=H{P + Qi) 
= P + Q, + Q^ 



(111.28) 



and generally 



(III.29) 



The above relations mean that the mapping H^{P) generates an additional subspace Q 



The sequence {P,H{P), ■ ■ ■ ,H'^{P)} is called the Krylov subspaces [13|. It may be clear 
that the subspaces {P, Qi, ■ ■ ■ , Qm} determine a unique block-tridiagonal form of H. In this 
sense the subspaces introduced in the present approach are essentially the same as those of 
Krylov. However, the basis states of each subspace Qk are ambiguous. Determination of the 
basis states depends on the purpose, that is, what problem we want to solve after the block 
tridiagonalization of the Hamiltonian. We show, in the later sections, that the basis states 
introduced in the present study are useful for the formulation of the effective-interaction 
theory. 



B. Expression of the Q box in terms of the uj operator 



Here we define two operators 



e{E)=Q{E-H)Q, 
1 



X{E) 



e{E) 



QHP 



1 



:QHP. 



E-QHQ 

In terms of x(-E), the Q box in Eq. (III.16]) is expressed as 



(IH.30) 



(III.31) 



Q{E) = PHP + PHQx{E) 



(IH.32) 



and the solution oj in Eq. flII.14l) to the decoupling equation flII.6P is given by 



UJ 



Y,x{Ek)\<Pk){'Pk 



(III.33) 
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where and (0^1 have been defined in Eqs. OII.lip and (III.12p . Consequently calculating 
Q{E) reduces to that of x{E). When the Q space is decomposed as in Eq. (IIII.17p . also xi^) 
is as 

X{E) = xi{E) + X2{E) + ■■■ + Xn{E) + ■ ■ ■ , (111.34) 

where 

Xn{E) = Q„ x{E)P. (III.35) 

Coupled equations for {xn{,E)} 

Qi e{E){xi{E) + X2{E)} = Q^HP, (III.36) 

e{E){xi{E) + X2{E) + xsiE)} = 0, (111.37) 

Qne{E){Xn-l{E)+Xn{E)+Xn+l{E)} = 0, (III.38) 

are derived from Eq. flIII.31D using Eqs.(ITlL7!), flIII.24D and flIII.25D . Since the Q box is 
expressed as 

Q{E) = PHP + PHQ^Xi{E) (III.39) 
by using Eq. (lIII.7p . calculating the Q box reduces to that of xi{E). 

C. Expansion in terms of continued fraction 



We show that the Q box is expanded by a continued fraction 40] of small-dimensional 



matrices by solving Eqs. (]III.36P - (1111.38^ . We assume Xm{E) = for m > 2 , then we have 

XiiE) = -^^QiHP (III.40) 



from Eq. ( jIII.3^ . where 

ei{E) = Q^{E - H)Q^. (III.41) 

Hereafter we use the notation 

e^{E) = Qm{E - H)Q^. (III.42) 
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The solution f lIII.40|) gives the Q box in the first approximation as 



Q^'\E) = PHP + PHQi-^QiHP. (III.43) 



Next we have 



X2{E) = -^Q2HQ,xi{E) (III.44) 

from Eq. (lIII.37p by assuming Xm{.E) = for m > 3. Substituing this into Eq.( IIII.36|l leads 

to 

Xi{E) = ^— Q,HP, (111.45) 

ei{E)-Q,HQ2——Q2HQ, 
e2{E) 

and then the Q box is given as 

Q(2)(^) = PHP + PHQi J QiHP (III.46) 

e,{E) - Q,HQ2—^Q2HQi 
e2{E) 

in the second approximation. Repeating similar manipulations, we finally have a general 
form 

Q{E) = PHP + PHQi ^ QiHP (III.47) 

ei — Hi2 J H21 

62 — H23 J Hi2 

63 — H34 Hi3 

64 

with Cm = em{E) and 

Hi, = QiHQj. (III.48) 

Here we consider a case in which the Q space for a system of interest is well described by 
finite number of subspaces. We denote the maximum of n by in Eq.f llll.lTl) . We introduce 



{en{E)} given through a descending recurrence relation starting from n = as 

en-i{E) = en-i{E) - Hn-l,n^Jr^Hn,n-l, (III.49) 

where we define 

eN{E) = Qn{E - H)Qn. (III.50) 
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H 



eXB 



H 



FIG. 3: Diagrammatical expression of Q{E) in terms of continued fraction. The thick hne denotes 
the renormalized propagator (ei(£'))^^. 



From Eq. (IIII.49|) we have a sequence eN^i{E), eM-2{.E), ■ ■ ■ , and ei{E). Then the Q box is 
expressed as 

Q{E) = PHP + PHQi,J—QiHP. (III.51) 

ei[E) 

Diagrammatical expression of Q{E) is shown in FigJS] It is a remarkable fact that the 
above result for the Q box indicates the existence of the renormalized inverse propagator 
ei{E) such that the Q box can be represented by a sum of only two terms, namely, the 
unperturbed part and the second-order term. 

If the dimension of the Q space is finite, the number of the subspaces {Qm} is also finite 
and the Q box given in Eq.f lIII.5T]) is exact. On the other hand, if the dimension of the 



Q space is infinite, the number of the subspaces {Qm} is, in general, infinite. For this 
case we introduce a truncation of the Q space. We consider a finite-dimensional subspace 
Qi + Q2 + ■ ■ ■ + Qn, where the subspaces {Qm, I < m < N} lead to a block-tridiagonal 
form of H as in Eq. (1III.26l) . The operator ei{E) that is determined through the recurrence 
relation in Eq.( ]III.4"9|) starting with n = is a function of and we write it as 'e{^\E). If 
e[^\E) converges as N tends to infinity, we can write the Q box as 

Q{E) = PHP + PHQ,::^^—Q^HP, (III.52) 

where 

= lim ■^,^\e). (III.53) 

A''— >-oo 
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D. Expansion in terms of renormalized vertices and propagators 



We here consider a method of calculation by an ascending recurrence relation for {xn{E)} 
and derive another solution for the Q box. By using Eqs. flIII.30p and flIII.48p . the coupled 
equations Eqs. (1111.360 - (]III.38p for the operators {xn{E)} are written as 

e,{E)xi{E) = H,o + HuX2{E), (III.54) 
e2{E)x2{E) = H2iXi{E) + H^^XsiE), (III.55) 

en{E)Xn{E) = H^,n-lXn^l{E) + Hn,n+lXn+l{E) , (III.56) 

with 

iJio = QiHP. (III.57) 

Equations (1III.56P is a linear relation of three operators Xn-iiXm and Xn+i) which can be 
cast into those of two operators as follows: First we rewrite (]III.54p as 

Xx{E) = a,{E) + ME)X2{E) (111.58) 

with 

a,{E) = -^H,,, (III.59) 

ME) = ^^12- (111.60) 



By substituting this into Eq. (jlll.SSp . X2{.E) is expressed as linear with xz{.E), 



where 



X2{E) = a2{E) + l32{E)xz{E), (III.61) 



a2{E) = ^ H21 I Hio 

/32(E) = J H23 

62 (-E) — H21 T— --^12 
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In general, the operators an{E) and Pn{E) obey the following ascending recurrence relations, 

1 



en{E) — Hn,n-lPn-l{E) 
f^n{E) = i—z — --^Hn,n+1, 



Hn,n-ian-l{E)^ 



and give a linear relation 



^n{.E) — Hn,n~ll3n^liE) 



Xn{E) = an{E) + Pn{E)Xn+l{E). 



(III.64) 
(III.65) 



(III.66) 



Equations ( 1111.641) and ( 1III.65I1 determine {a;„(-E'), (3n{E),n = 1, 2, ■ ■ ■ } with the initial values 
ai{E) and (5i{E) in Eqs. (]III.59[) and ( 1III.60I) . respectively. We finally have a solution for 

Xi(^) as 



Xi(^) = + Pi{E)a2{E) + ■■■ + Pi{E)P2{E) ■ ■ ■ /3„_i(E)a„(E) + 

CXI fe — 1 

= Y.[\[p„myk{E). 

k=l m=l 

Consequently the Q box is given by 



(III.67) 



Q{E)=PHP + PHQ^ 



oo k—1 



Y,{l[^rn{E)}ak{E) 



k=l m=l 



(111.68) 



In order to rewrite {(3n{E)}, and the Q box in terms of {ei{E)} and {Hij}, 

we introduce another inverse propagator em{E) defined through the following recurrence 
relation 



with the initial value 



(^miE) 6m{E) Hm,m—1_ / T-i\Hm—l,m 

em-i[E) 



ei(E) =ei(E) 

= Qi{E-H)Q^. 



(III.69) 



(III.70) 



We note that em{E) in Eq. ( 1III.6"9l) obeys an ascending recurrence relation, which differs from 
em{E) in Eq. (1III.49[) . In terms of {em{E)}, the operators {an{E)} and {f3n{E)} are written 
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as 



ai{E) = ^^//lo, (111.71) 

ME) = ^H,,^H,o, (III.72) 
e2{E) ei{E) 



Oin{E) = — r— r-f^n,n-lZ 7iqV-f^n-l,ra-2 " " " -^21 _ / -^10, (111.73) 

en{E) en-i{E) ei{E) 



ME) = ^^12, (111.74) 
HE) = ^^23, (111.75) 



PniE) = ^^Hn,n+l, (111.76) 



Then the Q box in Eq. 0III.68|) is expressed exphcitly as 



1 111 

Q{E) = PHP + Hoi _ , . Hio + Hoi _ Hi2 _ , . H2i _ Hio-\ 

ei{E) ei{E) e2{E) ei{E) 

+Hqi———Hi2 ■ ■ ■ „— — if„.„_i ■ ■ ■ H2i _ i/io + ■ • • . (III.77) 

ei[E) en{E) ei{E) 

A simpler expression of the Q box can be obtained by utihzing {Hk{E)} defined through 

= Hk-i{E)^^—Hk-i,k (III.78) 
ek^E) 

with the initial value 

Hi{E) = PHQi. (III.79) 
The Hk{E) interconnecting the P and Qk spaces is a dxdk matrix. The Q{E) in Eq. (11I1.77p 



is further reduced to 

QiE) = PHP + H,{E)^H\iE) + ■ ■ ■ + H^[E)^^h\{E) + ■ • • 

= PHP + Y,Hk{E)—-H[{E). (III.80) 
tit ^k[E) 
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This expression can be interpreted as that the Q box is given by a sum up to second 
order in the usual perturbation theory as schematically depicted in FigJH in terms of the 
renormalized inverse propagators {ek{E)} and the renormalized vertices {Hk{E)}. Equation 



0(£) = 



^ + Z o 



t = 1 



FIG. 4: Diagrammatical expression of Q{E) in the form of an expansion with the renormalized 
vertices and propagators. The thick line denotes the propagator (efc(i?))~^. 



(IIII.SOP shows clearly that there exist the renormalized inverse propagators {ek{E)} and the 
renormalized vertices {Hj^^E)} such that the Q box can be represented by a second-order- 
perturbation form which is the lowest-order interaction terms. If the original Hamiltonian H 
can be well described in a finite dimensional space, the Q box in Eq. (1III.80I) can be written 
by a sum of finite number of terms. 



IV. RECURSIVE SOLUTION FOR THE x{E) OPERATOR 

We here discuss how to calculate the operator {xn{,E)} in Eq.( lIII.35l) which are neces- 
sary for obtaining a true eigenstate, namely, |$fc) with the eigenvalue E = Ek- The basic 
equations for determining {xn{E)} have been given in subsections III B and III D. In the 
similarity-transformation theory for the effective interaction, the relationship between 1$^) 
and the model-space eigenstate 10^) is 

= \<l)k)+co\(pk). (IV.81) 
Using Eq. (]III.33P for u in terms of x(-^fc), l^fc) is also expressed as 

\<^k) = \(pk) + x{EkMk). (IV.82) 
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Therefore, if we want to obtain we have to solve x{Ek). We decompose x(-^fc) into 

{Xn{Ek)} as in Eq. flITL34D . The sequence {xi{Ek), X2{Ek), ■■■} obeys Egs- fHTLSiD -fHTrSGD. 
From Eq. (lIII.56p the following recurrence relation is obtained for {xn{Ek)}', 



Xn+liEk) — Kn-\-l^n{(^niEk)XniEk) H^^n—lXn—l [Ek)} {n > 2), 



(IV.83) 



where Kn,n+i is defined as 



n+1,71 



It is easy to see, using Eqs . (1111.231) and 0111.48^ for 



(IV.84) 



(IV.85) 



from which Eq. flIV.83P is derived. For the calculation of {xn{Ek)} with n > 3, Xi{Ek) and 
X2{.Ek) are necessary as initial values. In this stage we suppose that the Q{Ek) is given 
beforehand and use Eq. (1111.3^ to obtain 



with 



Xi{Ek) = K,oiQiEk)-PHP) 



(IV.86) 



(IV.87) 



where {pk\ and Iq^}^) are given in Eqs.f lIII.2|) and ( 1III.3p . respectively. In a similar manner, 
the operator X2{E) is solved, using Eq. (IIII.54|) . as 



X2iE) = K2i-{e,iE)xiiE)-H,,o} 



(IV.88) 



with 



(IV.89) 



where (g'^^''] and \q^^^) are given in Eqs. (]III.10|) and (lIII.lip . respectively. Substituting 
Xi{E) and X2{E), the sequence xsiE), X4:{E), ■ ■ ■ are obtained from the recurrence relation 
in Eq. fHv:83l) . 
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The eigenstate 1$^) of H with the eigenvalue Ek is finally given by 

\^k) = \<Pk) + Y,^n{Ek)\<Pk). (IV.90) 

n=l 

The usual normalization in the effective-interaction theory is {(f)k\4>k') = ^k,k'- Therefore, the 
normalized true eigenstate denoted by \^k) is given by 

\^k) = ^\^k), (IV.91) 



k 



where the normalization factor Nu is 



Nk = Jl + yiM^n{Ek)xn{Ek)\<pk). (IV.92) 



n=l 



V. THE Z-BOX METHOD AND EFFECTIVE HAMILTONIAN 



The Z box has been defined in the previous paper ^] as 



Z{E) = ^ [Q{E) - EQ,{E)] (V.l) 
1 - Qi{E) 



with 

dQiE) 



Qi{E) 



dE 



in order to overcome some defects that inevitably accompany the Q-box approach. The Z 
box has the following properties: 

(i) The operator 

d 

H,^ = Y,Z{Eu)\(t>k){4>k\, (V.3) 

k=l 

which is obtained by replacing Q{E) in Eg. ( lII.lTl) with Z{E), can be an effective 
Hamiltonian if {E^, = 1, 2, ■ ■ ■ , ci} are the eigenvalues of H. Therefore, Z{Ek) satis- 
fies the self-consistent equation 

Z{Ek)\<Pk) = Eu\<Pk)- (V.4) 
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(ii) The derivative of Z{E) is given by 

dZ{E) 2 



dE l-Qi{E) 



Q^{E)[Z{E) - EP] (V.5) 



with 



4w ^'''«^> 



2! dE^ 



Then 



dZ{E) 



dE 



E=E,, 



\4>k) = (V.7) 



holds for the eigenvalue E^ and the corresponding eigenstate 10^) of Hcs- 



(iii) For the eigenvalue Eg of QHQ determined by Eg. (111.20]) . Z{eq) satisfies the self- 
consistent equation 

Z{e,)\ii,)=sM. (V.8) 
Here we note that belongs to the P space. 

(iv) Contrary to Eq.f lV.7p . 

dZ{E) 



\n) = 21/^,) (V.9) 

E=ea 



dE 

holds for the derivative of Z{E) at E = Eg. 

These properties lead to conclusions: Z{E) is finite and differentiable even at £^ = e^, a 
pole of Q{E). Although E = Eg is also a solution of the self-consistent equation for Z(E), 
they can be easily discriminated from true eigenvalues {E^} of H with the aid of their 
derivatives, in Eqs. (]V.7I) and (]V.9I) . The Z-hox method has been applied recently to a 



2l| 



realistic calculation of the effective interaction by Coraggio et al. 

In order to calculate the Z box we need the first and second derivatives of the Q box. 
These derivatives can be calculated analytically and are derived in Appendices A and B 
corresponding to two expressions of the Q box given in Subsections III C and III D, respec- 
tively. 
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VI. MODEL CALCULATION 



A. Graphical method for eigenvalues of H 

We shall solve the eigenvalue problem for the Hamiltonian H in the framework of the 
Z-hox theory. We note that the Z box is a dimensional operator acting in the P space and, 
therefore, it has d eigenvalues. We have assumed that the operator Z(E) for an arbitrary 
energy variable E has d different eigenvalues. In the present calculation we do not discuss 
the case that Z{E) has some degenerate eigenvalues. The eigenvalues of Z{E) are functions 
of E. We write the eigenvalue equation for Z(E) as 

ZiE)\a) = Fk{E)\a), k = l,2,---,d. (VI.IO) 

The above eigenvalue equation defines d functions {Fk{E),k = 1,2, ■ ■ ■ ,d}. We label 
{Fk{E)} in order of energies as Fi{E) < F2{E) < ■■■ < Fd{E). From Eq. (1V.4I) we see 
that the solutions for the eigenvalues of H can be obtained by solving 

Fk{E) = E. (VI.ll) 

As shown in the previous Section, Eq.( ]V.4l) has two kinds of solutions, namely, E = Ei 



and E = Ej, where Ei and ej are the eigenvalues of H and QHQ, respectively. We distinguish 
the eigenvalues {ej} from {E^} according to the condition that the energy derivative dZ /dE 
takes different values for E = Ei and E = Ej. We define functions {Fl^{E), k = 1,2, ■ ■ ■ ,d} 
as 

FkiE) = (Ck ^ a), (VI.12) 



dE 

where \Ck) is the eigenstate given in Eq.f lVI.lOj) . The functions take the values 



= for E = Ei, (VI.13) 

Fl^{E) = 2 for E = Ej. (VI. 14) 

From the above properties of we see that the eigenvalues {Ei} of H can be obtained 

by calculating the solutions satisfying Eqs.l jVI.TTl) and (IVI.13I) simultaneously. A simple 
expression of the equation to be solved may be written as 

9.{E) = [^^^^^^]\{F^iE)r 

= 0, (VI. 15) 
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where Fq is a parameter chosen suitably such that the two terms on the right-hand side take 
values of the same order of magnitude. 

The solutions to Eq. (IVI.15|) can be obtained by a graphical method. We define a function 
fk{E) as 

where A is a small number. The function fk{E) has properties 



1 



hm^ fk{E) = — (VI.17) 



and 



for the eigenvalues Ei of H and ej of QHQ, which may be obvious from Eqs. ( IVI.TTTl -f lVI.TSl) . 
If the parameter A is taken to be small enough, the function fk{E) behaves like a resonance 
at E = Ei. By drawing the graph of {fk{E), k = 1,2, ■ ■ ■ ,d} and finding resonance positions, 
we obtain eigenvalues of H. 

B. Numerical calculation 

In order to obtain some assessments of the present approach we study a model problem. 
We start with a model Hamiltonian H of which matrix elements are given by 

{i\H\j) = (m + Pi^)5i, + 7X,, (VI.19) 

with 

- 1, (VI.20) 



Xij 



2|^v^(z+j) 



where [X] is Gauss's symbol which means the integer part of a real number X. A set of 
{xij} are recognized as pseudo random numbers satisfying 

-1 < Xij < 1. (VI.21) 

The a, /3 and 7 are the parameters chosen suitably. The total dimension of H is taken to be 
Nfi = 100. As for the P space we choose a two-dimensional space {d = 2) spanned by the 
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two states which have the lowest and second lowest diagonal energies of H. We here do not 
consider a case that some of the eigenvalues {y^*"^} in Eqs. fllll.lOl) and fllll.2ip become zero, 
because {xij} are pseudo random numbers and H does not have any definite symmetry. 
Therefore, the subspace {Qk, k = 1,2, ■■■ ,Nq} are all (i-dimensional and the number of 
the subspace {Qk} is given by Ng = {N^ — 2)/2 = 49. 

We first calculate the Q box and its energy derivatives Qi{E) and Q2{E) according to 
the continued-fraction method and the renormalized vertex method formulated in Sections 
111 C and 111 D, respectively. We have confirmed numerically that the calculations on these 
two methods agree to each other. With Q{E), Qi{E) and Q2{,E) we calculate the Z box 
and its energy derivative dZ{E)/dE according to Eqs. (lV.ip and (IV. Sp . 

We next calculate the functions Fk{E) and given in Eqs. flVl.lOl) and flVl.121) . 

respectively. We finally obtain the functions {fk{E), k = 1,2, ■ ■ ■ ,d} and draw graphs of 
these functions. Since the dimensions of the P space is taken to he d = 2, we have two 
graphs of fi{E) and f2{E). These graphs are shown in FigJSl From these figures we can 
specify the position of the eigenvalues of H as resonance positions. From FigJSl we can 
estimate four eigenvalues of H on the interval [0, 10]. 

5 
4 



o 
O 




FIG. 5: Resonance-like behavior of the functions fi{E) and f2{E) in the case of a = 1.2, (3 
0.2, 7 = 1.4, A = 10^2 and Fq = 1.0. 



The accurate solution, namely Ei, can be obtained in the following way: We suppose that 
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the solution Ei lies on the interv al [a , h] and there are no other solutions on this interval. The 



parabolic-interpolation method 4l| is applied here. If the difference \E — Ei\ is sufficiently 
small, the approximate form of gk{E) in Eq.f lVI.l^ becomes a parabolic function written as 

9,iE) = i±i^oW)i!(i^ _ E,r. (VI.22) 

Therefore we approximate gk{E) to be a parabolic function and solve the energy Ei to give 
the minimum of gk{E). We note here that the parabolic function A{x — a)^ passing through 
two points {a,gk{a)) and {b,gk{b)) takes the minimum at the point a given by 

a = ^vgg+^vgg, (VL23) 



where we have assumed a < a < b. We utilize this fact to solve Eg. (1VI.15I1 . 



The calculation procedure employed in this numerical calculation is as follows: 

(i) Determine an interval [a, b] on which only one solution Ei exists. 

(ii) Divide [a, b] into equal intervals and define five points {Ei, E2, E3, E^^, E^) as 

Ej, = a + {k- 1)AE, l<k<5 (VI.24) 

with AE = (6-a)/4. 

(iii) Consider all the intervals [Ei, Ej] by selecting Ei and Ej among {Ei, E2, ■ ■ ■ , E^} and 
calculate 

^rUEi) + ^/UE~) 

(iv) There are ten combinations of the energies {Eij}. Arrange {Eij} in order of energy 
and write them as Ui < M2 < ■ • ■ < Uiq. 

(v) Calculate the values {gk{un)-, n = 1, 2, ■ ■ ■ , 10} and find the minimum gk{um) as shown 
in Figjni We determine a new interval [a, b] according to 

a = Um~i, b = Um+i if 2 < m < 9, 
a = a, b = U2 if m = 1, 

a = ug, b = b if m = 10. (VI.26) 
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(vi) Repeat the procedure until the convergence, \gk{um)\ < S, is attained for an appropri- 
ate small number 6. 



fk{E) 




a U1U2 ■ ■ ■ Um-lU 



Um+1 ■■ - Uio b 



» E 



FIG. 6: Determination of a new interval. If fk{um) is the minimum value among {fkiun), n = 
1, 2, • • • , 10}, then the new interval is given by [a, b] = [um-i,Um+i]- 

In Table [T] we show the results for the lowest two eigenvalues of H calculated by the above 
mentioned parabolic-interpolation method. The convergence is markedly fast. With three 
times of the changes of the interval [a, b], convergence is reached with accuracy better than 
10 decimal places. 

As has been shown in Eq. (IIII.8"0|) . the Q box is given by a sum over the number k. 
In this model calculation the maximum number of k is equal to Nq = 49. Introducing a 
number -ft"max, we consider a truncation as k < -ft'max in the calculation of the Q box. It 
would be interesting to examine the dependence of the calculated eigenvalues of H on -ft'max- 
The results are shown in Figs |7] and [H It is clear that, as -ft'max approaches to Nq=49, the 
eigenvalues converge to the exact values. These results suggest a possibility of introducing 
a truncation in the series expansion for the Q box. 



VII. CONCLUDING REMARKS 



We have proposed a new approach to the effective interaction and/or Hamiltonian acting 
within a model space P. In the present stage of the effective-interaction theory one of the 
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TABLE I: Correct digits of the lowest two eigenvalues of H calculated by the parabolic-interpolation 
method. The parameters a, /?, 7 are taken as the same as in FiglSj 

Initial intervals are taken to be [a,b] = [0.0, 1.0] and [2.5,3.5] for Ei and E2, respectively. 



Ei No. of repeats Calculated value 

El 1 0.365- • • 

2 0.365550- • • 

3 0.36555015199 
E2 1 2.999- • • 

2 2.9994240- • • 

3 2.9994240873 



0.55 

0.50 
1 0-45 

a 

V 
M 

s 

0.40 



"""'0 10 20 30 40 50 

Block Dimension ATmax 

FIG. 7: Convergence of Ei as a function of i^max- 

central problems has been how to calculate accurately the Q box which has been used as 
a building block of the formulation. The main concern of the present study has been to 
derive a new method of calculating the Q box as accurately as possible even if the original 
Hamiltonian H is given in a huge-dimensional space. 

The formulation consists of two steps: First one is to transform a given Hamiltonian 
H to a block-tridiagonal form by dividing the complementary space Q of the P space into 
subspaces {Qk, k = 1,2, ■ ■ ■} with tractable dimensions. If the subspaces are chosen suitably 
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FIG. 8: Convergence of E2 as a function of 



the Hamiltonian is transformed to a block-tridiagonal form. With the Hamiltonian thus 
transformed, the next step is to derive coupled equations for determining the Q box. By 
solving these coupled equations we have proved that the Q box can be represented in two 
ways : (a) The first one is that the Q box is expanded into a form of continued fraction in 
terms of the submatrices which are the elements of the block-tridiagonalized Hamiltonian. 
It has been proved that if a quantum system can be well described by a Hamiltonian given 
in a finite dimensional space, the continued fraction can be reduced to only one term with a 
renormalized propagator which can be calculated by using a descending recurrence relation. 
In this case the number of the subspaces {Qk, k = 1,2,---} should be finite, (b) The 
other solution is obtained by using ascending recurrence relations for solving the coupled 
equations. The resultant Q box can be shown to be given by only two terms such as PHP 
and a sum of second-order terms with respect to renormalized vertices and propagators. This 
reduction of the Q box has clarified that there exists a method of determining renormalized 
vertices and propagators such that the Q box, or equivalently the effective Hamiltonian of 
the Bloch-Horowitz or the Feshbach form, can be given by a sum of terms up to second 
order. 

Given the Q box, we have applied the Z-hox method for solving the eigenvalue problem 
of a Hamiltonian H. We have introduced functions of energy variable E as {fk{E), k = 
1, 2, ■ ■ ■ ,d} such that fk{E) has a pole at E = Ei which is one of the eigenvalues of H. Here 
the number d is the dimension of the model space. In this approach the eigenvalues of H 
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can be given by the pole positions of the functions {fk{E), k = 1,2, ■ ■ ■ , d}. This approach 
enables us to solve the eigenvalue equation of in a graphical way. 

In order to assess the present method we have made a test calculation by introducing 
a 100x100 model Hamiltonian. We have performed the calculation of the Q box by em- 
ploying two methods, namely, the continued-fraction expansion and the expansion with the 
renormalized vertices and propagators. We have confirmed that both the two methods have 
reproduced the exact eigenvalues of the original Hamiltonian H. 

The present study is based essentially on the algebraic approach to the effective Hamilto- 
nian. If we want to calculate the effective interactions acting among valence particles outside 
the core, it is necessary to represent the Q box in terms of linked diagrams. Therefore, a 
remaining problem is how to relate the renormalized vertices and propagators introduced in 
the present study with the linked-diagram approach. 
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Appendix A: Derivatives of the Q box in Eq. (jIII.51|) 

The first and second derivatives of the Q box are given, respectively, by 

where ei{E),ki{E) and li{E) are given through the following recurrence relations: We 
consider the energy derivative of en{E) in Eq.( lIII.4"9l) and write as 



k^E) = (A.3) 
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Noting a relation 



1 denjE) 1 



= (A.4) 

we can derive 

1-1 

kn{E) = Qn + Hn,n+1- -f^kn+l{E)- —Hn+l,n, (A.5) 

en+i{E) en+i{E) 
where we have used the energy derivative of en{E) in Eq. (lIII.4ip 

^^Q.. (A.a) 

The Hn^n+i and i^„+i,n are defined in Eq. flIII.48l) . For the maximum number of n, denoted 
by A^, k^i^E) is given by 

~ deN^E) 
kn{E) = 

= Qn, (A.7) 



which is derived from Eq. piI.50P for ei^{E). Starting with kfyf{E), the recurrence relation 



determines a sequence kiy{E), kiy_i{E), ■ ■ ■ , ki{E). 
We write the second derivative of e„(i?) as 

d^er,{E) 



ln{E) 



dE^ 
dkn{E) 



dE 

From Eq. ( lA.Sp for {/c„(i?)} a recurrence formula for {/„(£')} can be derived as 

ln{E) = — 2Hnn+l— TTW^n+liE)— ——^n+l{E)— -——Hn+l,n 

en+i[E) en+i[E) e„+i(E) 



(A.8) 



+ Hn,n+1- j^ln+l{E)- -—Hn+l,n- (A.9) 



For the maximum number n = N the In{E) is given, from Eqs. (lA.7p and (lA.Sp . by 



In{E) = 0. (A.IO) 

The recurrence formula Eq. (lA.9p determines a sequence In{E),In-i{E), ■ ■ ■ ,li{E). Substi- 
tuting the operators ei{E),ki{E) and li{E) into Eqs. llA.ll) and (\A.2\} the first and second 
derivatives of the Q box can be calculated. 
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Here it should be noted that the first and second derivatives of the Q box can be expressed 
by using only small- dimensional matrices. The ei(i?), ki{E) and li{E) are the operators on 
the subspace Qi, and they are represented by di^di matrices, where di is the dimension of 
the Qi space. The operator PHQi is a mapping between the P and Qi spaces and has a 
d^di matrix representation, where d is the dimensions of the P space. 

Appendix B: Derivatives of the Q box in Eq. (IIII.80|) 

We derive the first and second derivatives of the Q box with respect to energy variable 
E as 

= f;{(:ffl(i5;)A.(i5;):^i(E) + h.c.) +:ff,(E)AUi5;):ffI(i5;)}, (b.i) 

fc=i 

^^^^ = Y.{2%{E)\,{E)H'l{E) + {hI{E)\,{E)hI{E) + h.c.) 



fc=i 



with 



+ 2{H[{E)\mHl{E) + h.c.) + H,{E)\1{E)hI{E)} (B.2) 



%iE) = (B.3) 
HliE) = (B.4) 
XkiE) ^ (B.5) 



KiE) 



d\k[E) 
dE 

1 dek{E) 1 

ek{E) dE ek{Ey 

dE^ 

1 dckiE) 1 dckiE) 1 1 d%{E) 1 



(B.6) 



(B.7) 



ek{E) dE ek{E) dE ek{E) ek{E) dE^ ek{E)' 
These expressions indicate that the calculation of the derivatives of Q{E) is reduced to 
that of {Hk{E)}, {Xk{E)} and their derivatives; {Hk{E)} are given through the recurrence 
relation Eq.f llll.TS!) and accordingly its derivatives are 

H[iE) = {H',_,iE)\k-iiE) + Hk-iiE)\',_,iE)}Hk-,,k, (B.8) 
HliE) = {h'1_,{E)\,^,{E) + 2%_,{E)\',_,{E) + H,^,{E)\U{E)}H,^r,,, (B.9) 
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and {\k{E)} is given by 



from Eqs. flIII.69p and fIB.Sp . and accordingly its derivatives are 
■^U-^) = ~ '^ki^) + ^k{E)Hk,k-i\'k-i{E)Hk-i^k\{E) , 



(B.IO) 



(B.ll) 



K{E) = - {X'j,{E)Xk{E) + h.c.) + {X',{E)Hk,k-iK-i{E)Hk-i,kXk{E) + h.c.} 

+ Xk{E)Hk^k-iXl_i{E)Hk-i^kXkiE). (B.12) 



Their initial values are given in Eq. (IIII.79p for {Hk{E)}, H^{E) = H^{E) = and 

1 



Xi{E) 



ei{E) 



{E - QiHQj 



X[{E) = -{E - Q.HQi] 
X'liE) = 2{E - QiHQ, 



\-3 



(B.13) 
(B.14) 
(B.15) 



The way of calculating the derivatives of the Q box is summarized as follows: First, 
{Xk{E)} is calculated by Eq. llB.lOp . then their derivatives by Eqs.f lB.lT]) and f lB.12p . Next 
{H'k{E)] and {H"k{E)} by Eqs.(E8]) and f lR9|l . finally we obtain the derivatives of the Q 
box. Here Xk{E), A'^(i?), and X'l{E) are d^xdh matrices, while Hk{E)^ H'^{E), and h'1{E) 
are dxdk matrices. 
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